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ABSTRACT. The evolution of two point space cor- 
relation function of QSOs {S,QQ{r, z)) is analyzed in the 
framework of the theory of the large scale structure for- 
mation. For given cosmological models the agreement 
between theoretical predictions and observational data 
on ^QQ{r,z) is achieved by determination of scales of 
structures R{z), in which could be formed quasars at 
different redshifts z. It is shown that quasars could be 
formed in the nucleus of massive galaxies or groups of 
galaxies owing to merges at their central region. It 
depends on an amplitude of the power spectrum of 
the initial matter density fluctuations at galaxy-cluster 
scales. 
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1. Introduction 

In the problem of the large-scale structure forma- 
tion the QSOs present a special interest because they 
cover almost all epoch of evolution of the Universe af- 
ter cosmological recombination. The most physically 
motivated explanation of quasar phenomena is the ac- 
cretion of gas and stars into the massive black holes in 
the center of galaxies. It would take place at the early 
stage of galaxy evolution or later as a result of collision 
or merging of galaxies in the groups or clusters. There- 
fore, the theoretical precalculation of spatial two-point 
correlation functions of QSOs [S.QQi''', z)) has an im- 
portant meaning for testing the cosmological scenarios 
of the formation of the large-scale structure of the Uni- 
verse. 

The analysis of quasar surveys (lovino et al. 1991; 
Andreani, Cristiani 1992; Mo, Fang 1993; Romberg, 
Kravtsov 1995; Fang 1998; La Franca et al. 1998; 
Cristiani et al. 1998 and references therein) indicates 
on the clustering of QSOs higher than galaxies and 
lower than rich clusters. The observable quasar cor- 



relation function can be approximated by the expres- 
sion ^QQ{r,z) = (r/r^)-^(l + z)"" with 7 « 1.8 and 
w 6 — 10/i~^Mpc. As it was shown by lovino et 
al. 1991, Mo, Fang 1993 and Romberg, Rravtsov 1995 
the amplitude of ^gqir, z) decreases with increasing of 
redshift z {a > 0). But they give the different rate 
(values of a) of evolution. Andreani, Cristiani 1992 
did not reveal an essential evolution of CgQ(?', z). 

Contrary, taking into account the last data on QSO 
space distribution La Franca et al. 1998 and Cristiani 
et al. 1998 show that amplitude of ^QQ(r, z) increases 
with z. According to their results 

eQ%{l5h-\z) = (ro«/r)i-«(l + z)-(i-2+^), (1) 

where the "fitting-parameter" e = —2.5 ± 1.0. Here 
^QQ{15h~^Mpc, z) is integrated correlation function 
defined as 

r 

£.QQ{riz) ^ ^ j x^^QQ{x,z)dx. (2) 


The goal of this paper is theoretical interpretation of 
such evolution of ^qq(15/i~^A/pc, z) within the frame- 
work of cosmological models with the given initial 
power spectra P{k) and existing assumptions about 
physical nature of quasars. 

2. Quasars correlation function 

We assume that the QSOs are manifestation of short- 
term active processes, which occur in nucleuses of mas- 
sive galaxies at the early stage of their evolution or 
later one owing the merging of galaxies in their groups 
or clusters (Haehnelt 1993; Haehnelt, Rees 1993; Loeb, 
Rasio 1994; Rontorovich 1996; Efstathiou, Rees 1998). 
According to the most elaborated scenario of the large 
scale structure formation the galaxies, their groups and 
clusters are formed in the peaks of the random Gaus- 
sian field of the matter density linear perturbations 
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owing to gravitational instability. So, the amplitude 
of relative density fluctuations S = = v ■ (j{R) [v 
is the height of peaks, (j{R) is the root mean square 
of density fluctuation at corresponding scale R) are 
normally distributed. The quasar stage in an evolu- 
tion of such structures begins at some moment of time 
t (or redshift z), when massive black hole at center 
of galaxies is formed or close galaxies in most dense 
region of groups or clusters are merged. We can ne- 
glect an interval of time At, during which the quasar 
mechanism is formed, because it is essentially smaller 
than cosmological time t (At ^ t) as it was shown 
by Efstathiou, Rees 1988, Haehneh 1993, Haehneh, 
Rees 1993, Loeb, Rasio 1994 and Kontorovich 1996. It 
means that the moment of beginning of quasar stage 
is close to the moment of appearence of the singular- 
ity in the Tolman solution for the collapse of dust-like 
cloud. In the initial random field of the density fluc- 
tuations the peaks of corresponding scale and ampli- 
tude result into such collapsing clouds. The values of 
5 for fluctuations which coUaps at redshift z is equal 
to 5{z) — 5c' {z + 1). In a case of spherical - symmet- 
ric fluctuations - 5c — 1.69, nonspherical - 5c = 1-33 
(Efstathiou, Rees 1988, Nusser, Silk 1993). 

In the framework of the theory of the random Gaus- 
sian fields (Bardeen et al. 1986) the two-point cor- 
relation function of peaks in which the QSOs will be 
formed is following: 



^QQ{f, z) = {^\liQQ{r.z) + yjis5{r,z) 



(3) 



where r - distance between them in comoving coordi- 
nates, £,qq{t, z) is the statistical component of correla- 
tion function of peaks in which the QSOs on given z are 
formed, £^ss{r,z) is the correlation function of matter 
density fluctuations, it is the dynamical component. 

For model with given initial power spectrum of den- 
sity fluctuations P{k) the correlation function ^i5i5(r, z) 
is equal 



oo 

6.(^^) = ^-/fc^-^ 



P{k,R) sinikrl 



where P{k,R) = P{k) ■ W^{kR), W{kR) = 
exp (— 5 • fc^ • R^) is the Gaussian window function for 
the scale R, which corresponds to scale of fluctuations. 
The height of peaks of density fluctuations, in which 
the QSOs will be formed on fixed z, is defined as: 



uiz) 



5c-{l + z) 
aiR) 



Then ^55 (r, z) can be presented as follows: 

2 



^ss{r,z) 



v{z) 



(5) 



(6) 



where x{r) = £,ss{r)/^ss{0), £,ss{r) = ^ss{r, z = 0). 

The expression for £^qq (r, z) was obtained by 
Novosyadlyj, Chornij 1998 by the method like Kaiser 
1984 one for correlation function of rich clusters but 
taking into account that QSO phenomenon is short- 
term stage after collapse (or merging) in center of cor- 
responding proto-object. It has the following view: 



^QQir,z) 



1 



• exp 



i^^(z) • x{r) 
1 + x{r) 



1. (7) 



Substituting expression (||),(@) in expression (||), we 
have complete expression for correlation function of 
quasars £,QQ{r, z) which we analyze here. 

3. Evolution of correlation function of quasars 

At the beginning let's consider the case when the 
QSOs are formed in structures of the same fixed scale 
R, which does not depend on cosmological time t. For 
this we will analyze the derivative of ^qq (r, z) with 
respect to redshift z, ^^^^j^^, at any fixed r, and 
find conditions which define its sign. So long as the 
matter density correlation function in the range where 
< x{r) < <7^{R) (that corresponds for most models 
r < 30 — 40/i~^Mpc) is monotonous function of r its 
derivative with respect to redshift z does not depend 
on value of r. Therefore without losing of generality 
we can analyze this derivative analytically in the range 
x{r) <C 1 where the expression (|^) becomes simpler: 



(8) 



The QSO correlation function (||) in approach x{r) <C 1 
will look as follows: 



iQQ{r,z) ^ v{z) + 



v{z) 



x{r). 



(9) 



It means that 



dz 



= when 



v{z*) = V<5c. (10) 

Taking into account (H) we shall obtain: 

a(i?) = V^-(z, + 1). (11) 

It means, that in the case when QSOs are formed 
in the peaks of the same scale R at different z 
their correlation function £,QQ{r,z) monotonously in- 
(^dioQ^ > 0) if v(z) > and decreases 

^ ^ ^) ^^^^ V^- The interesting case 

is when condition ( [lO[ ) occurs at > 0. In this 
case the amplitude of £,QQ(r, z) decreases in the range 
< z < z*, reaches its minimum at z = z* and in- 
creases at z > z*. 



creases 

( diQQ{r,z) 
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Figure 1: Izoline of minima of amplitudes of QSO 
correlation function S.QQ{r,z) in the plane <j{R) — z. 
By filled circles it is shown the moments at which the 
evolution of QSO correlation function changes the sign 
for for standard CDM model, tilted ones and stan- 
dard H+CDM (30%Hot+60%CDM+10%baryon) for 
the case when quasars are formed in the peaks of galaxy 
scale R ~ 0.35/1" ^ Mpc. 



In particular, in the case when the amplitude of 
^QQ(r, z) decreases with increasing of z, or in the case 
of absence of essential evolution, that was indicated 
by lovino et al. 1991, Andreani, Cristiani 1992, Mo, 
Fang 1993, Komberg, Kravtsov 1995, it may mean that 
QSOs are formed in low peaks of a random Gaussian 
field of matter density fluctuations, ^{z) < y/S^. All 
these cases are shown in Fig.l. The solid line divides 
the plane into the ranges with decreasing (left) and 
increasing (right) the amplitude of QSOs when z in- 
creases. The moments at which the evolution changes 
the sign for standard CDM model, tilted ones and stan- 
dard H-hCDM (30%Hot+60%CDM-|-10%baryon) one 
are shown by filled circles. 

The values of for CDM models with slopes of 
their power spectra n — 1, 0.8, 0.7 are 3.16, 1.06, 
0.45 accordingly For H-fCDM model = 0.21. 
Thus, in the interval < z < 3.16 the amplitude of 
^QQ(r, z) monotonously decreases only in CDM model 
with n = 1 that agree with the data by lovino et 
al. 1991; Andreani, Cristiani 1992; Mo, Fang 1993; 
Komberg, Kravtsov 1995. However, as it is shown 
by Novosyadlyj, Chornij 1998 owing to small height 
of peaks the amplitude of the correlation function of 
quasars £,QQ{r, z) in this case becomes less than galac- 
tic one already at z ^ 1, that does not match men- 
tioned above observational data. It may mean that 
quasar phenomena would be connected with structures 
of larger scale than galaxy one. 



Therefore, let's analyze the case when scale of ob- 
jects in which QSOs appear is different at different 
redshifts and determine them using the last data on 
the evolution of QSOs correlation function S,QQ{r,z) 
given by La Franca et al. 1998 and Cristiani et al. 
1998. It can be realised in the following way: a) we 
set some values for cosmological parameters that give 
us possibility to calculate the power spectrum of den- 
sity perturbations using the analytical approximation 
of H-I-CDM power spectra in four dimensional param- 
eter space (Novosyadlyj et al. 1999), b) we deter- 
mine the R{z) using the fQg(15ft.~^Mpc, z) given by 
as input observable data at different z. Thereby we 
can determine the typical scales of structures R{z), in 
which the quasars at different z were formed. The scale 
i?(z) can specify conditions and mechanism of forma- 
tion of QSOs: when i?^0.5 — lft.~^ Mpc it suggests 
that main part of QSOs appeared in the nuclei of mas- 
sive galaxies which are not dynamically connected with 
structures of larger scales, when R> 1 ^ 2h~^ Mpc the 
formation of quasars in groups of galaxies due to pro- 
cesses of merges of galaxies is more probable. 

The calculations were carried out for redshifts in 
the range 0.5 < z < 3.5 and matter dominated flat 
Universe with following values of the cosmological pa- 
rameters: rim = 1 (matter density in units of critical 
one), content of baryons ilf, = 0.096 followed from nu- 
cleosynthesis constraint (Tytler et al. 1996), content 
of the neutrino hot dark matter Q^, — 0.1 (one sort 
of massive neutrino), content of the cold dark matter 
^CDM = ^m — ^b — ^}^, dimensionless Hubble constant 
h = Ho/100{km/ sec/Mpc) = 0.5, and different slopes 
of the primordial power spectrum - n — 0.8, 1, 1.2. The 
initial power spectra normalized to 4-year COBE data 
(Bennet et al. 1996, Bunn & White 1997) are shown 
in Fig. 2. For each flxed value of redshift z we calcu- 
lated the integrated correlation function of QSOs (2) 
using the formulas (3, 4, 7), equated it to correspond- 
ing observational one (1) for r = 15/i~^Mpc and solved 
the obtained equation relative to R{z). The results for 
three models are shown in Fig. 3. 

As one can see, in the models with low amplitude of 
the power spectrum at small scales {n = 0.8) the QSOs 
are formed in the peaks of galaxy scale R = 0.5^ IMpc 
{M ^ IO^Mq). And contrary, if the amplitude of the 
power spectrum is higher at small scales (n = 1 — 1.2) 
then the preferable range of scale of structures in which 
QSOs appeared are R = 2 ^ 3Mpc, that corresponds 
to galaxy groups and clusters (M ^ 10^^ 10^^ Mq). 
It can be caused by merging of galaxies in their most 
dense regions at the center. We do not discuss here the 
drastic decreasing of such scales at small redshifts for 
the model with ti = 1.2 because the analysis presented 
here has rather the qualitative character than quanti- 
tative. For the more certain quantitative predictions of 
scales of structures connected with QSOs phenomena 
the more realistic physical models of their formation 
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Figure 2: Power spectra of density fluctuations P{k) 
at z = (4- year COBE normalization) for cosmological 
models with parameters flm = 1, = 0.096, fli, = 
0.1, ilcDM = 0.804, h = 0.5, and different slopes of 
primordial power spectrum n = 0.8, 1, 1.2. 




of dynamical one. It can be provided by choosing 
of corresponding scale of structures R{z) in which 
the QSOs were formed at different z . 

• The QSOs can be formed as in nuclei of single 
massive galaxies as in consequence of merging of 
galaxies in their groups or clusters. Which mech- 
anism is statistically dominant it depends on an 
amplitude of the initial power spectrum at small 
scale. 
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4. Conclusions 

Based on the results presented here we conclude: 

• The growth of the amplitude of the two-point 
space correlation function of QSOs i^QQ{r,z)) 
with increasing of z, which was revealed by La 
Franca et al. 1998; Cristiani et al. 1998 can be 
explained within the framework of the theory of 
formation of the large-scale structure of the Uni- 
verse from the peaks of the random Gaussian field 
of matter density fluctuations. The matter of ex- 
planation is that growth of statistical component 
of ^QQ (r, z) should be faster than the reduction 
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